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Abstract 

The results of an experimental calibration of the NASA 
Glenn Research Center 16 in. Mass-Flow Plug (MFP) are 
presented and compared to a previously obtained calibration 
of a 15 in. Mass-Flow Plug. An ASME low-beta, long-radius 
nozzle was used as the calibration reference. The discharge 
coefficient for the ASME nozzle was obtained by numerically 
simulating the flow through the nozzle from the WIND-US 
code. The results showed agreement between the 15 and 16 in. 
MFPs for area ratios (MFP to pipe area ratio) greater than 0.6 
but deviate at area ratios below this value for reasons that are 
not fully understood. A general uncertainty analysis was also 
performed and indicates that large uncertainties in the 
calibration are present for low MFP area ratios. 

Introduction 

Large-scale model supersonic aircraft inlets are routinely 
tested in NASA Glenn Research Center’s (GRC) larger wind 
tunnels such as the 8- by 6-fl and 10- by 10-ft Supersonic 
Wind Tunnels (SWT). Two key objectives of such tests are to 
ensure that the air supplied by an inlet matches the mass-flow 
demand of an engine and that the inlet operates in a stable 
condition over a wide range of mass-flows. Thus, accurate 
measurement of inlet mass-flow is of primary importance. Due 
to its relatively simple mechanical design, the mass flow 
through an inlet is often controlled and measured by a 
calibrated conical mass-flow plug (MFP). For larger models, 
GRC uses a 16 in. diameter MFP which was fabricated by 
modifying a pre-existing 15 in. MFP. Calibration of the 16 in. 
MFP plug was performed in-house using an ASME low-beta, 
long-radius nozzle as the reference measurement. The 
procedure and resulting calibration, as well as a comparison to 
the calibration of the previous 15 in. MFP are reported herein. 


A Area (in’) 

AR MFP area to cold-pipe area ratio (Eq. (20)) 

C D Discharge coefficient 
D Diameter (in) 

L i ASME nozzle throat length (in, Figure 3) 

M Mach number 
MFP Mass Flow Plug 
R Radius (in) 

Rj ASME nozzle ellipse major radius (in. Figure 3) 

R 2 ASME nozzle ellipse minor radius (in. Figure 3) 

R air Gas constant for air (53.3525 ftTbf/(lb-R) 

Re,; ASME nozzle Reynolds number 
Re,;^ Scaled ASME nozzle Reynolds number 
Re</, s = Re,; x lOx 10~ 6 

T Temperature (R) or throat tap location (in, Figure 3) 

U Velocity (ft/sec) 
d ASME nozzle diameter (in) 
g c Proportionality constant (32.17405 lb-ft/(lbf-s 2 ) 
p Pressure (psi) 

s Slant height of the frustum of a cone (in. Figure 5) 

w mass flow rate (lb/s) 

x,r, 9 Cylindrical coordinates 

x,y,z Cartesian coordinates 

Xp| ug Measured position of MFP (in) 

x'p| U g Effective position of MFP (in. Figure 5 and Eq. (23)) 

P Ratio of ASME nozzle-to-approach pipe diameter 

5V, Uncertainty of measure and X t 

4> MFP cone half-angle (deg, Figure 5) 

y Ratio of specific heats 

p Density (lb/ft 3 ) 
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Subscripts 

a Pertaining to the average component 
amb Pertaining to ambient conditions 
i Pertaining to ideal conditions 
n Pertaining to the normal component 
noz Pertaining to the ASME nozzle 
pipe Pertaining to the cold-pipe 
plug Pertaining to the MFP 
pot Pertaining to potential flow conditions 
t Pertaining to total conditions 
v Pertaining to the vertical component 
5* Pertaining to boundary-layer blockage (displacement 
thickness) 

Superscripts 

* Pertaining to choked flow conditions 


Pam b (^) 

T amb<^> 



Calibration Procedure 

The calibration of both the caged 15 and 16 in. mass-flow 
plugs were carried out with the experimental setup shown in 
Figure 1. The four primary components of the calibration rig 
are the American Society of Mechanical Engineers (ASME) 
inlet flow nozzle assembly, an adapter spool piece, the cold- 
pipe assembly, and the mass-flow plug assembly. The entire 
calibration assembly was mounted vertically on the floor of 
GRC’s 10- by 10-ft Supersonic Wind Tunnel and protruded 
through a sealed opening in the ceiling. The ASME nozzle 
was fed by atmospheric ambient pressure and temperature air 
while the wind tunnel exhausters were used as a vacuum 
source to establish a flow through the assembly. A photo of 
the ASME nozzle on the top of the 10- by 10-ft. SWT is 
shown in Figure 2 and a schematic of the geometry is shown 
in Figure 3. A photo of the 16 in. MFP attached to the cold- 
pipe is shown in Figure 4 and a schematic of the plug is shown 
in Figure 5. 

The spool piece is an adapter to match the ASME nozzle 
exit diameter to the cold-pipe diameter under test. The cold- 
pipe is simply a 96 in. long constant-area pipe whose primary 
purpose is to allow for a flow-settling length between the inlet 
exit and the mass-flow plug. Near the forward end is a 
spherical-nosed centerbody supported by three equally spaced 
struts. The centerbody nose is intended to simulate the fan 
spinner on a turbojet engine. Attached to the forward end of 
the cold-pipe is a spool piece that houses eight compressor 
face rakes. The rakes lie in the plane of the tip of the 
centerbody nose. The inside edge at the exit plane of the cold- 
pipe is sharp to fix the position of the mass-flow plug choke 
point. 



Figure 2. — ASME nozzle assembly. 
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A picture of the 16-in. mass-flow plug attached to the cold- 
pipe is shown in Figure 4. The cone is actuated by a hydraulic 
cylinder and concentricity of the cone with the cold-pipe is 
maintained by wedge bearings on three of the six knife-edged 
support struts (cage). A position transducer (potentiometer) 
can be seen on the top of the assembly in Figure 4. 

The ultimate goal of the calibration is to determine the 
discharge coefficient for the 16 in. MFP over an operating 
range typical of the supersonic inlet testing. The discharge 
coefficient for the mass-flow plug is defined as: 

r _ w noz 

t-'Zyplug — * ' ‘ ) 

Wplug,; 

where w noz is the actual mass flow through the cold-pipe and 
Wpiug,,- is the ideal inviscid choked flow through the mass- 
flow plug. For the purpose of calibrating the mass-flow plug, 
the actual mass flow was measured with an ASME Long- 
Radius Flow Nozzle. A description of this nozzle is as 
follows. 

ASME Nozzle 

Geometry: An ASME Long-Radius Flow Nozzle with 
throat taps was used to provide a reference mass-flow measure 
for calibrating both the 15 and 16 in. mass-flow plugs. The 


geometry of the nozzle is shown in Figure 3. This design 
conforms to the Low-p Nozzle with Throat Taps illustrated in 
Fig. 11-111-14 of Reference 2 with the following two 
exceptions. First, no approach pipe exists before the nozzle 
(D= co), hence (3=c//Z)=0. And second, the nozzle exit flow does 
not exhaust into a sudden expansion but rather into a shallow 
angle conical diffuser. The nozzle throat-to-conical diffuser 
and conical diffuser-to-constant area end section are both 
transitioned by 27.96 in. radii. 1 



Figure 4. — Caged 16 in. Mass-Flow Plug. 


Note that this diffuser geometry differs from the description given by Iek 
(Ref. 1) where the diffuser is described as being a “cubic contour.” 
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Citing the desire to maintain fully turbulent flow in the 
nozzle, the calibration tests were performed both with and 
without a boundary-layer trip in the inlet region of the nozzle. 
Details of the trip were not supplied by lek (Ref. 1) but 
inspection of the ASME nozzle indicated the trip location and 
width are as shown in Figure 3. Most of the tests were 
performed with an 80 grit trip, but limited testing was also 
performed with a 36 grit trip with essentially the same results. 

Flow Equations: The flow through the ASME nozzle is 
calculated by the following equation: 

w noz ~ C D,noz ' W noz,i (2) 


where 


T = 


T 

1 t,noz 

1 + 1— 1 ^--M 2 
1 ~ 2 1V± noz 


(7) 


into the expression for w 0 and rearranging yields: 


w„. 


Y 'gc 


Pt,noz ' ^ noz ' ^noz 


R a 




1 + 1 — --M 2 

1 ~ _ 1V± noz 


-(Y+l) 
- (y-i) 


( 8 ) 


Substituting the values for y, g c and R :m and replacing p Lnoz and 
Ft, noz with Pa m b and T amh , respectively, gives the final form: 


^noz.i ~ P noz * U noz ' ^noz 


(3) 


and where p noz and U„ oz are evaluated at the throat of the 
ASME nozzle. Substitution of the ideal gas law, Mach number 
definition, and isentropic and adiabatic relations for pressure 
and temperature, respectively: 


n _ 144 'Pno, 

(4) 

K-air ' * noz 


M n oz ' ’ Sc * -^air ' ^ noz 

(5) 

Pt,noz 

y 

(6) 



y-1 

2 


• M f W z 



w noZ j - 0.91884 • Pamb Aoz ■ (l + 0.2 -M„2 0Z )— 3 (9) 

V F am b 

where M noz is the Mach number of the flow at the throat of the 
nozzle and is given by the isentropic relation: 



f 


r i ) ^ 



r \ 


5- 


Pamb 

(3.5) 



V Pnoz y 



K 


/ 


( 10 ) 


and p noz is the static pressure measured at the throat of the 
nozzle. All of the variables in Equations (9) and (10) are 
directly measured or calculated from directly measured 
quantities. The only remaining parameter required to calculate 
the nozzle flow rate is the nozzle discharge coefficient ( C Dnoz ). 
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ASME Nozzle Discharge Coefficient 

Historical Correlation: The discharge coefficient for a 
nozzle with throat taps is a function of the potential flow non- 
uniformity and the boundary-layer displacement thickness 
resulting from viscous effects. As discussed by Rivas and 
Shapiro (Ref. 3), if the length of the cylindrical sections before 
and after the throat taps are on the order of or greater than l A 
throat diameter, then the effects of potential flow non- 
uniformity are negligibly small. 

In the absence of a calibration, ASME recommends reading 
the discharge coefficient for a low-(t long-radius nozzle with 
throat taps from Fig. II-III-19 of Reference 2. The source of 
this curve is not referenced but it appears to be the “Calculated 
Curve” in Fig. 22 of Cotton and Westcott (Ref. 4). No 
uncertainty is given for using this curve but based on similar 
analysis, it is likely to be on the order of ±1.5 to 2.0 percent. 
Much of the data used to generate this curve was based on 
incompressible flow data. In 1973, Smith and Matz (Ref. 5) 
investigated an 8-in. low-P, long-radius nozzle under 
incompressible and compressible conditions and reported a 
Mach number dependence. 

Two-Code Computational Approach: For the March 1996 
calibration of the 15 in. caged mass-flow plug (Ref. 1), due to 
the uncertain effect of the (S parameter and the nozzle exit 
diffuser on the discharge coefficient, the ASME correlation 
was not used. Instead, the theoretical method first performed 
by Smith and Matz (Ref. 5) and later by Lahti and Hamed 
(Ref. 6) was used to establish a discharge coefficient. With 
this method, the discharge coefficient is assumed to depend on 
two components, an inviscid part ( C DiPOt ) that accounts for the 
radial pressure variation in the throat, and a viscous part 
(Cad*) that accounts for the effective boundary-layer 
blockage. Computer programs were used to calculate each 
component and the discharge coefficient of the nozzle was 
determined from the product of the components: 

Cd,iioz, 2C — Co, pot ‘Cafi* (1 1) 

This approach will hereafter be referred to as the “Two- 
Code” (2C) nozzle calibration. In Lahti and Flamed’s (Ref. 6) 
study, the potential solution was obtained with the Stream 
Tube Curvature (STC) computer code (Ref. 7), chosen 
because of their many years of experience with it and because 
it is driven by global mass conservation as part of its solution 
procedure making it a natural choice for mass-flow calibration 
studies. The boundary-layer correction was obtained using 
Harris and Blanchard’s code (Refs. 8 and 9). This code was 
chosen based on its wide use and thorough documentation. 
Lahti and Hamed (Ref. 6) report that comparisons with 
experimentally determined discharge coefficients indicate that 
an accuracy of ±0.25 percent can be achieved with their 
method. 


For the calibration of the GRC ASME nozzle, lek (Ref. 1) 
followed a similar approach but used different computer 
codes. For the potential solution, Stockman and Farrell’s code 
(Ref. 10) was used and for the viscous solution, Herring’s PL2 
code (Ref. 11) was used. With these codes, the theoretical 
discharge coefficient for the ASME nozzle was determined at 
six discrete throat Reynolds numbers. 


TABLE 1.— DISCHARGE COEFFICIENT 
VERSUS REYNOLDS NUMBER 


Re rfj =Rc rf xlO' 6 

^ D,noz 

1.6828 

0.98553 

3.0386 

0.98667 

3.7659 

0.98728 

4.0302 

0.98728 

4.3062 

0.98755 

4.6056 

0.98796 


For data reduction during the mass-flow plug calibrations, 
the data in Table 1 were interpolated with a third-order 
polynomial function: 2 

C D ,noz,2C = 0.982982 ± 0.00208 1573 Re^ - 

0.0004063900 Re^± 0.00004062582 Re^ s ' ' 

and extrapolated by setting points outside the calibration range 
to the end points of the calibration. 

RA NS Computational Approach: More recently, a different 
computational approach was taken for obtaining the discharge 
coefficient of the ASME nozzle as well as the 16 in. MFP. 
Rather than use a potential code coupled with a boundary- 
layer calculation, WIND-US, which is a general purpose 3-D 
Reynolds-averaged Navier-Stokes (RANS) code, was used for 
the numerical studies (Ref. 12). Versions 3 Alpha and 2 of 
WIND-US calculated the CFD derived discharge coefficients. 
WIND-US is a mature, multi-zone, structured-grid, 
compressible flow solver offering a variety of turbulence 
models. All simulations were run steady-state using the 
Menter Shear Stress Transport (SST) turbulence model 
(Ref. 13). The SST turbulence model was chosen for all of the 
simulations because previous studies (Refs. 14 and 15) have 
shown that this particular model performs well for jet flows. 

An axisymmetric structured grid, composed of roughly 
119,000 grid points and divided into 17 zones, was used to 
model the test rig assembly shown in Figure 1. Grids were 
created for 19 different plug positions, ranging from x p i ug = 2.4 
to 16.4 in. with x p i ug defined as shown in Figure 5. Zones 
upstream of the experimental boundary-layer trip location (see 
Figure 3) assumed laminar flow while zones downstream of the 


2 

Reference 1 indicates a fourth-order polynomial was used but the actual data 
reduction program uses the third-order fit shown. 
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trip location assumed turbulent flow. An example grid is shown 
in Figure 6 while Figure 7 illustrates the solid computer model. 

The simulations used a constant Courant-Friedrichs-Lewy 
(CFL) number of 0.5 rather than a constant time step. To 
reduce computational costs and convergence time, parallel 
processing along with three levels of grid sequencing were 
used to obtain converged solutions. First, a coarse grid 
sequence included only every fourth point in each direction. 


Second, a medium grid sequence was used that utilized every 
other point in each direction. Third, a fine sequence used 
every grid point. Solution convergence was achieved when the 
computational equation residuals were reduced by several 
orders of magnitude and the mass flow rates in the bellmouth 
entrance and near the cold pipe exit each varied by less than 
0.09 percent per iteration. Convergence was typically obtained 
in a little over 3 hr. 


100 - 


c 

> 


50 - 


-100 


-50 



0 50 

X (in) 

Figure 6. — Sample primary grid layout, laminar zones are in blue, turbulent zones are in red. 



Figure 7. — The 90° sector of the axisymmetric solid model. 
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Initially, the W1ND-US code was run with the as-tested 
ambient inlet conditions (p am t= 14.3 psia), but then to better 
resolve Reynolds number effects, additional cases were run at 
inlet pressures of 5.0 and 50.0 psia. For all pressures, the 
ambient total temperature was set to 520.0 °R with an initial 
inflow Mach number set to 0.001. Whereas for the 5.0 psia 
case, the back-pressure at the MFP exit was set at 1.0 psia, for 
the two higher pressure cases, the back-pressure was set at the 
as-tested pressure of 2.0 psia. 

After the CFD analyses were completed, the data was post- 
processed to acquire the discharge coefficients for the 
bellmouth nozzle. First, the mass flow rate through the cold 
pipe was calculated from the average mass flow rate between 
the downstream tangential edge of the radial nose of the 
conical plug to the cold pipe exit: 


N 

!■ 

i = 1 


pipc.x.; 


N 


(13) 


where 

Wpipe,*=jjp -U x -r-dr-dQ (14) 

(() r 


This was compared to the average mass flow rate at the 
throat of the bellmouth to check for continuity. The solution 
was converged when the average mass flow rates in the 


bellmouth throat and just upstream of the cold pipe exit varied 
by less than 0.09 percent of each other. The discharge 
coefficient could then be calculated from: 


C D , 


K pipe 


(15) 


where w nozJ is given by Equation (9). 

The data from the W1ND-US computations are plotted in 
Figure 8. For reference, the ASME curve, the data from Smith 
and Matz, and the numerical data from lek are also shown. For 
the WIND-US data, only points where the throat Mach 
number was between 0.2 and 0.8 are plotted. Due to the use of 
a boundary-layer trip on the nozzle, the distribution does not 
display the transitional behavior of Smith and Matz (Ref. 5) 
data and indeed seems to agree well with the turbulent portion 
of their “Critical” curve (Rep„ ( , : > 2.0xl0 6 ). Another 
observation of the current data is the saw-tooth nature of the 
combined distributions indicating a small Mach number 
dependence. To illustrate this dependence more clearly, data at 
constant Mach numbers of 0.3, 0.5 and 0.7 were extracted and 
are plotted in Figure 9. Although the dependence on Mach 
number is small for the current data, the trend is opposite that 
of Smith and Matz (Ref. 5), but since their low Mach number 
data was in the transition region, nothing firm can be drawn 
from this. Comparing the present WIND-US data to the 
previous distribution established for the same nozzle 
configuration by lek (Ref. 1) shows that the present data 
trends about 0.2 percent higher than Iek’s. 



Figure 8. — ASME nozzle discharge coefficient. 


N ASA/TM— 20 14-218128 


7 


1.005 


1.000 

0.995 

J- 0.990 

U 

0.985 

0.980 

0.975 

























* 





X 












- 






















X 

V 

V 



















^ o 

c 

D a 

JP 



























ASME (Ref. 2) O lek (Ref. 1) 

S&M, M~0.4 (Ref .5) — ■— WIND, M-0.30 

S&M, M~0.7 (Ref.5) — WIND, M=0.50 

- • - S&M, Critical (Ref.5) — •— WIND, M=0.70 

, i i ■ 1 — ' ' 1 1 — 












0.1 1.0 10.0 


Re d (X10 6 ) 


Figure 9. — ASME nozzle discharge coefficient, constant Mach number. 


In order to apply easily the present data to the 16 in. MFP 
calibration, a correlation was developed that fits the data and 
accounts for the small Mach number dependence. This was 
accomplished by first curve fitting the M=0.5 data which 
resulted in the following expression: 


Flow Equations: The discharge coefficient for the mass- 
flow plug is assumed to be a function only of the plug position 
and defined as the ratio of the actual mass-flow to the ideal 
theoretical mass-flow (see Eq. (1)). The ideal mass-flow is 
calculated from the following equation: 


Cd,M= 0.5 


= 0.993003- 


6.73896 

sfRe d 


(16) 


An empirical correction to this nominal discharge 
coefficient was established such that the Mach number 
dependant discharge coefficient is defined by: 


where 


Wplug,! 


" Ppipe ' ^pipe ' ^R ' ' ^plug 


ar = 

^pipe 


(19) 


( 20 ) 


Cd,hoz,WU - Cd,M= 0.5 ' f( M ) (17) 

where 

/(M)= 0.999514 • (l - 0.2M2 oz )“ 000863415 (18) 

Since it is unlikely that the current ASME nozzle would ever 
be operated at anything but ambient pressure conditions, the 
coefficients in Equation (18) were determined from the 
14.3 psia data only. A comparison of Equation (17) and 
Equation (16) with the numerical data is shown in Figure 
10(a) and Figure 10(b), respectively. 


and where A p iug is the choked-flow area between the cold-pipe 
exit plane edge and the conical plug. With reference to Figure 5, 
the definition of the choked area (/l p | llg ) of the MFP is based on 
the area of the frustum of a right circular cone where the slant 
height (s a ) is defined as the average of the vertical distance from 
the cold-pipe edge to the MFP and the slant height of a cone 
whose slant height is normal to the MFP surface. 

dplug = tl ' $a (^pipe + R a j (21) 

where (see Figure 5) 


16 in. Mass-Flow Plug 

Geometry: With reference to Figure 5, the plug itself is a 
4> = 16.5° half-angle cone with a 2.18 in. nose radius. The 
downstream end of the cone is 7? p i ug =15.800 in. in diameter 
followed by a 0.610 in. wide cylindrical shoulder. 


s n — -*piug • sin(<)>) 

S v — “^piug ' tan((|>) 
s a — (s n + S v )/2.0 

R a = -Kpipe - s a ' COs((f) - areCOs(.S'„ /s a )) 


( 22 ) 
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(b) Fit of all data. 


Figure 10. — Comparison of ASME discharge coefficient correlation with numerical data. 


and 


-'■plug ~ -'-plug + (^pipe - -Kplug )/ tan (<l>) ( 23 ) 

where x ' p i ug is the distance from the end of the cold-pipe to the 
theoretical intersection line of the plug conical surface and the 
cylinder of the cold-pipe internal diameter (see Figure 5). The 
choice of s„ to approximate the sonic line (curved line in 


Figure 5) is arbitrary but must be used consistently for the 
calibration to be valid. 

The area ratio can be expressed in terms of the Mach 
number in the pipe by the following isentropic relation: 


^plug 
A ■ 

^pipe 


y+1 

y + l) 2 (y-l) 


M 


pipe 



y-i 

2 


■M 2 - 

pipe 


~(Y+1) 
2(Y— 1) 


(24) 
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Substituting Equation (24) and Equations (4) to (7) (with the 
subscript “ noz ” replaced with “pipe”) into Equation (19) 
yields: 


W/.plug 



— (Y+l) 

Y + lV(y-l) 
2 , 


(25) 


The total pressure in the cold-pipe (p t: p i pe ) is calculated from 
Equation (6) (with the subscript “«oz” replaced with “pipe”). 
The static pressure in the cold-pipe (p p j pe ) is directly measured 
with static pressure taps and the Mach number in the cold-pipe 
{M p ipe) is calculated from Equation (24) which assumes that 
the flow at the mass-flow plug exit is choked. 

Instrumentation: With reference to Figure 1, the test setup 
was instrumented as follows. The ambient air pressure and 
temperature were each recorded with four transducers which 
were then averaged to obtain p amb and T amb . There were a total 
of six throat taps on the ASME nozzle and eight static taps at 
the mid-point of the cold pipe. These were also averaged to 
obtain p noz and /; plpe . The temperatures were measured with 
Type-E thermocouples and the pressures were measured with 
ESP 8400R 0 to 15 psia modules. The data was recorded on 
the GRC ESCORT data acquisition system. 

Calibration Runs: For the caged 16 in. mass-flow plug 
calibration runs, there were six sweeps of the mass flow plug. 
The first three were performed with the boundary-layer trip 
installed on the ASME nozzle and the last three were 
performed with the trip removed. 


well the CFD models the flow physics. In particular, since the 
nozzle discharge coefficient is primarily a function of the 
boundary-layer growth, the turbulence model will have an 
influence on the result. Further, the difference between how 
the boundary layer transitions on the nozzle will also affect the 
discharge coefficient. In the CFD, the flow is assumed laminar 
until the leading edge of the boundary-layer trip and fully 
turbulent downstream. In the experiment, there is a finite 
width of grit material. These effects are not easily quantifiable, 
thus the uncertainty specified for the nozzle discharge 
coefficient is to be considered a “best guess” by the authors. 
Finally, although the ambient temperature is measured and 
used in the calculations, it does not contribute to the MFP 
discharge coefficient uncertainty. 


TABLE 2.— MEASURAND UNCERTAINTY 


i 

A,- 

8 Xt 

(15-/1 6-in. MFP) 

Units 

i 

■^•plug 

0.015 

in. 

2 

T 

1 amb 

2.5 

°R 

3 

Pamb 

0.0060/0.0045 

psi 

4 

Pnoz 

0.0049/0.0037 

psi 

5 

P pipe 

0.0042/0.0032 

psi 

6 

d 

0.0025 

in. 

7 

^pipe 

0.0075 

in. 

8 

Dpiug 

0.0050 

in. 

9 

<l> 

0.05 

deg 

10 

C D, noz 

0.001 

N/A 


Uncertainty Considerations 

A general uncertainty analysis was performed for both the 
15 and 16 in. MFP calibrations. The uncertainty in the MFP 
discharge coefficients was estimated by: 



where X, are the measurands which are summarized with their 
assumed uncertainty in Table 2. 

The measurands P amb , P, to - and / J plpe are actually calculated 
averages of four, six and eight individual pressure 
measurements, respectively, and the uncertainties in Table 2 
reflect the measurement redundancy. In addition, different 
pressure uncertainty levels were used for the 15 and 16 in. 
MFP calibrations to reflect an upgraded pressure measurement 
system for the 16 in. MFP. The measurement uncertainty for 
the individual modules was 0.012 and 0.009 psi for the 15 in. 
and 16 in. MFP, respectively (Ref. 16). The uncertainty in the 
nozzle discharge coefficient is comprised of two components. 
The first is related to the grid and solution convergence and is 
generally small (<0.003 percent). The second is related to how 


Results and Discussion 

Caged 15 in. MFP Calibration 

The 15 in. MFP was calibrated over a range from 
Xpiug = 0.25 to 17.0 in. with and without distortion screens 
installed in the cold-pipe. A correlation was established based 
on the data without distortion screens. The discharge 
coefficient was assumed to be only a function of the choked 
flow area ratio (Eq. (20)) and is expressed as a segmented 
third-order polynomial function where the coefficient values 
are given in Table 3: 

Q),piug,2C = c o + c i ' + ci ■ AR 2 + c 3 ■ AR^ (27) 


TABLE 3.— COEFFICIENT VALUES FOR 
CAGED 15 in. MASS-FLOW PLUG 



AR < 0.203943 

0.203943 <AR< 0.81 1515 

0.81 15 15 < AR 

Co 

1.0340732 

0.96163863 

1.1845804 

Cl 

-0.58958996 

0.0047836495 

-0.60737883 

C 2 

1.1527497 

-0.030427921 

0.42665471 

C 3 

0 

0.050854378 

0 
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This original calibration (Eq. (27)) was based on the two- 
code ASME nozzle calibration. To update the calibration to 
use the more recent W1ND-US based nozzle calibration, two 
approaches are possible. First, the data calibration data points 
could be recalculated with the W1ND-US based nozzle 
calculation and a new fit established, or second, an adjustment 
could be applied to Equation (27) by multiplying by the ratio 
of nozzle discharge coefficients: 

Q),plug,»X/ = Co.plug.2C ■ {C D ,tioz,WU /C D , noz ,2c) (28) 

The nozzle discharge coefficients based on both methods 
are plotted in Figure 11 as a function of the 15 in. MFP area 
ratio. Also plotted is the ratio of the coefficients which can be 
approximated by the following empirical relation: 

Q 

D ' noz ’ WU = 1.00 1 4446 + 0.0 1 7657 1 09 • AR 

C D,noz,2C 

- 0.00069852959/ AR - 0.0468 18586 ■ AR 2 (29) 

+ 2.3567944 x10" 5 /AR 2 +0.049185331- AR 3 

- 4.4263225 x 10“ 7 /AR 3 - 0.0182043 18 • AR 4 

The original and adjusted 15 in. MFP calibrations along 
with the data they were derived from are shown in Figure 12. 


The nozzle Mach number and Reynolds number are also 
plotted for reference. With regard to the original calibration, 
there are a couple of things to note about this distribution. At 
area ratios below approximately 0.2, the curve has a large 
negative slope and at the lowest area ratio, the discharge 
coefficient exceeds unity with large data scatter. At the largest 
area ratio, the slope and scatter also increases. During the test, 
the ratio of the cold-pipe base pressure to pipe total pressure is 
monitored to ensure the plug is operating under choked 
conditions. The two highest area ratio positions are 
approaching flow conditions that challenge the assumption of 
a well-defined choked flow. 

The adjusted 15 in. MFP discharge coefficient curve fit is 
re -plotted in Figure 13 with the uncertainty bounds from the 
general uncertainty analysis. Also plotted is the percent 
uncertainty based on the assumed measurand uncertainties 
(Table 2). 

These results indicate that the measurement uncertainty 
exceeds 0.5 percent for MFP area ratios below 0.5. Based on 
this calibration, the inlet should be sized to operate within the 
area ratio range 0.2<AR<0.8 (1.85 in.<jt plug <l 0.4 in.). 

The relative contribution from the individual measurands is 
plotted in Figure 14. Also plotted for reference is the absolute 
uncertainty. The large uncertainty at low area ratios is 
primarily due to the low Mach number in the nozzle (see 
Figure 12) which is calculated from Equation (10). 
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Figure 1 1 . — ASME Nozzle discharge coefficient as a function of 15 in. MFP area ratio. 
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Figure 12. — Original and adjusted 15 in. MFP discharge coefficient 
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Figure 13. — The 15 in. MFP discharge coefficient uncertainty. 


N ASA/TM— 20 14-218128 


12 


% Uncertainty Nozzle Mach or Re (xl.0E-07) Number 




Figure 14. — Measurand contribution to overall 15 in. MFP discharge coefficient uncertainty. 


Caged 16 in. MFP Calibration 

The results of the calibration of the 16 in. MFP for both 
tripped and un-tripped nozzle boundary layer are shown 
Figure 15. The nozzle Mach number and Reynolds number are 
also plotted for reference. At a given area ratio, due to its 
larger size, the 16 in. MFP has a larger flow rate than the 15 
in. MFP. At large area ratios, this results in choking of the 
nozzle for the 16 in. MFP. Thus data beyond an area ratio of 
0.7 will be excluded from further analysis as the MFP would 
be operating unchoked. 

Comparing the 16 in. MFP data with and without the 
boundary-layer trip shows agreement above an area ratio of 
about 0.5, but deviates significantly below this value. At an 
area ratio of 0.5, the nozzle Reynolds number is about 
3.0xl0 6 . From Figure 9, a Reynolds number of 3.0xl0 6 is 
consistent with the end-of-transition implied by the data of 
Smith and Matz (Ref. 5) and the ASME curve (Ref. 2). Smith 
and Matz data show a large increase in discharge coefficient 
below this transition value. Since a fully turbulent nozzle 
discharge coefficient was used in the 16 in. MFP calibration, 
which was artificially low for the case with sweeps with no 
boundary-layer trip, the resulting MFP calibration is under 
predicted. In fact, what the transitional nozzle discharge 
coefficient would have been could be estimated by correcting 
the no-trip data to agree with the tripped data. 

A comparison between the adjusted 15 in. MFP and the 
tripped 16 in. MFP calibration data is shown in Figure 16. The 
16 in. MFP calibration data above an area ratio of 0.7 has been 


removed (ASME nozzle choked) and the remainder of the data 
was fit with a second order polynomial: 

C Aplug = 0.943044 + 0.026915 ■ AR + 0.017376 ■ AR 2 (30) 

Also included in this plot are the uncertainty bounds from 
the general uncertainty analysis. The experimentally 
determined discharge coefficients for the 15 and 16 in. MFP 
tend to agree well above an area ratio of 0.6, but deviate below 
this ratio. In particular, the 16 in. MFP curve does not exhibit 
the abrupt increase in flow coefficient below an area ratio of 
0.2. Inasmuch as the calibration curves lay at the edges of the 
others uncertainty bounds, this suggests that the differences in 
calibration are not fully attributable to measurement 
uncertainties. For both MFPs, the method of calibration 
produces large uncertainties at low area ratios. This is 
primarily due to the low Mach number in the ASME nozzle. 
The use of a delta pressure transducer at low Mach numbers 
would mitigate the large uncertainties significantly. 

As previously mentioned, the numerical study with WIND- 
US included a numerically based determination of the 16 in. 
MFP discharge coefficient (see Figure 6 and Figure 7). This 
study is still in progress and will be reported in a future 
publication, but the preliminary results indicate a trend very 
similar to the experimental 16 in. MFP results but shifted up 
about 0.25 percent. Flowever, results for area ratios below 0.2 
have not yet been obtained so no comment can be made 
regarding the tail-up behavior of the 15 in. MFP. 
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Figure 15. — The 16 in. MFP discharge coefficient data. 



Concluding Remarks 

The results of a calibration procedure for the GRC 16 in. 
MFP have been presented and compared to a previous 
calibration of a 15 in. MFP. The 16 in. MFP was fabricated 
from the 15 in. MFP so all future testing will be with the 16 in. 
MFP. For the 16 in. MFP, a calibration has been established, but 
a difference in character between the 15 in. MFP and the 16 in. 
MFP for area ratios below 0.6 is observed. Unfortunately, this 
region of disagreement coincides with large uncertainty levels. 
The recommended path forward is to first complete the 


numerical calibration of the 16 in. MFP; in particular obtaining 
calibration data below an area ratio of 0.2. And second, the 
experimental calibration of the 16 in. MFP should be repeated 
with a focus on minimizing measurement uncertainties. 
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